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In ultrarelativistic heavy-ion colhsions, the Fourier decomposition of the relative azimuthal an- 
gle, A<^, distribution of particle pairs yields a large cos(3A0) component, extending out to large 
rapidity separations A77 > 1. This component captures a significant portion of the ridge and shoul- 
der structures in the At;!) distribution, which have been observed after contributions from elliptic 
flow are subtracted. An average finite triangularity due to event-by-event fluctuations in the initial 
matter distribution, followed by collective flow, naturally produces a cos(3A(/>) correlation. Using 
ideal and viscous hydrodynamics, and transport theory, we study the physics of triangular {vz) flow 
in comparison to elliptic (^2), quadrangular (va) and pentagonal (ws) flow. We make quantitative 
predictions for U3 at RHIC and LHC as a function of centrality and transverse momentum. Our 
results for the centrality dependence of W3 show a quantitative agreement with data extracted from 
previous correlation measurements by the STAR collaboration. This study supports previous re- 
sults on the importance of triangular flow in the understanding of ridge and shoulder structures. 
Triangular flow is found to be a sensitive probe of initial geometry fluctuations and viscosity. 



I. INTRODUCTION 

Correlations between particles produced in ultrarela- 
tivistic heavy-ion collisions have been thoroughly studied 
experimentally Correlation structures previously iden- 
tified in proton-proton collisions have been observed to 
be modified and patterns which are specific to nucleus- 
nucleus collisions have been revealed. The dominant fea- 
ture in two-particle correlations is elliptic flow, one of the 
early observations at RHIC [l|. Elliptic flow leads to a 
cos(2A0) term in the distribution of particle pairs with 
relative azimuthal angle Ac/). More recently, additional 
structures have been identified in azimuthal correlations 
after accounting for contributions from elliptic flow. 0t3| ■ 
An excess of correlated particles are observed in a narrow 
"ridge" near A0 = and the away side peak at A0 = tt 
is wider in comparison to proton-proton collisions. For 
central collisions and high transverse momentum triggers, 
the away side structure develops a dip at A^ = tt with 
two "shoulders" appearing. These ridge and shoulder 
structures persist for large values of the relative rapidity 
A?7, which means that they are produced at a very early 
times 3. 

It has been recently argued [01 that both the ridge and 
the shoulder are natural consequences of the triangular 
flow (W3) produced by a triangular fluctuation of the ini- 
tial distribution. The purpose of this paper is to carry 
out a systematic study of v^ using relativistic viscous hy- 
drodynamics, which is the standard model for ultrarela- 
tivistic heavy- ion collisions [10]. We also perform trans- 
port calculations [1 Ij , because they allow us to check the 
range of validity of viscous hydrodynamics, and also be- 
cause they provide further insight into the physics. Along 
with U3 , we also investigate W4 (quadrangular flow) and W5 
(pentagonal flow). In Sec. [Hi we recall why odd moments 
of the azimuthal distributions, such as V3, are relevant. 



In Sec. lIIIi we study the general properties of anisotropic 
flow induced by a harmonic deformation of the initial 
density profile using hydrodynamics and kinetic theory. 
In Sec. IIVI we present our predictions for v^ and U5 at 
RHIC and LHC. The contribution of quadrangular fluc- 
tuations to Vi is difficult to evaluate because Vi also has 
a large contribution from elliptic flow 12]: this will be 
studied in a forthcoming publication [13i] . 



II. CORRELATIONS FROM FLUCTUATIONS 

A fluid at freeze-out emits particles whose azimuthal 
distribution /(0) depends on the distribution of the fluid 
velocity [l^l- f[<P) can generally be written as a Fourier 
series 
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where w„ are the coefficients of anisotropic flow [14] which 
are real and positive, and i/'n is defined modulo 2-k /n (for 
fn 7^ 0). Equivalently, one can write 
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where angular brackets denote an average value over out- 
going particles. 

Generally, «„ is measured using the event-plane 
method [15|. However, two-particle correlation measure- 
ments are also sensitive to anisotropic flow. Consider a 
pair of particles with azimuthal angles 0i , (/'2 = 0i + A0. 
Assuming that the only correlation between the particles 
is due to the collective expansion, Eq. ([2]) gives 
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The left-hand side can be measured experimentally, and 
Vn can thus be extracted from Eq. ([3]) [16j. Experimen- 
tally, one averages over several events. Vn fluctuates 
from one event to the other, and the observable mea- 
sured through Eq. ([3]) is the average value of (w„)^. It 
can be shown that the event-plane method also measures 



transverse plane (x, y) is a deformed Gaussian at t = to 



the RMS, 



unless the "reaction plane resolution" is 



extremely good [17], |18| . 

Most fluid calculations of heavy-ion collisions are done 
with smooth initial profiles |19l423l |. These profiles are 
symmetric with respect to the reaction plane f/'i?, so that 
all V'n in Eq- © are equal to f/'fl (with this convention, 
all Vn are not necessarily positive). For symmetric colli- 
sions at midrapidity, smooth profiles are also symmetric 
under — >■ -|- tt, so that all odd harmonics wi, v^, etc. 
are identically zero. However, it has been shown that 
fluctuations in the positions of nucleons within the col- 
liding nuclei may lead to significant deviations from the 
smooth profiles event- by-event [2J, [25|. They result in 
lumpy initial conditions which have no particular sym- 
metry, and this lumpiness should be taken into account 
in fluid dynamical calculations .26ji29l|. More precisely, 
one should calculate the azimuthal distribution for each 
initial condition, then average over initial conditions. 

Initial geometry fluctuations are a priori important for 
all Vn, as anticipated in Ref. [30]. Their effect on flow 
measurements has already been considered for elliptic 
flow V2 [31J, Is^l and quadrangular flow va, [S^. Event- 
by-event elliptic flow fluctuations have been measured 
and found to be significantly large, consistent with the 
fluctuations in the nucleon positions [3J|. Directed flow, 
vi , is constrained by transverse momentum conservation 
which implies '^PtViipt) = and will not be considered 
here. In this paper, we study triangular flow vz [9|, and 
pentagonal flow v^, which arise solely due to initial ge- 
ometry fluctuations. 



III. 



FLOW FROM HARMONIC 
DEFORMATIONS 




FIG. 1: (Color online) Contour plots of the energy density 
(0) forn = 3 and £3 = 0.2. 

Elliptic flow is the response of the system to an ini- 
tial distribution with an elliptic shape in the transverse 
plane (cc, y) [S^]. In this article, we study the response to 
higher-order deformations. For sake of simplicity, we as- 
sume in this section that the initial energy profile in the 
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where we have introduced polar coordinates x — r cos 0, 
y = rsin0. In Eq. Q, n is a positive integer, e„ is the 
magnitude of the deformation, ij^n is a reference angle, 
and p is the transverse size. Convergence at infinity im- 
plies < e„ < 1. Fig. [Tj displays contour plots of the 
energy density for n = 3 and £3 — 0.2. The sign in front 
of En in Eq. ((4|) has been chosen such that tpn is the di- 
rection of the flat side of the polygon. For n = 2, it is the 
minor axis of the ellipse, which is the standard definition 
of the participant plane 25] 

For t > to, we assume that the system evolves accord- 
ing to the equations of hydrodynamics or to the Boltz- 
niann transport equation, until particles are emitted, and 
we compute the azimuthal distribution /(0) of outgoing 
particles. The initial profile @ is symmetric under the 
transformation — > -|- — , therefore /(</>) has the same 
symmetry. The only nonvanishing Fourier coefficients are 
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where u„ is real. As we shall see below, w„ is usually 
positive for £„ > 0, which means that anisotropic flow 
develops along the flat side of the polygon (see Fig. [T|) 

We now present quantitative results for w„, as defined 
by Eq. ([5|), using two models. The first model is rela- 
tivistic hydrodynamics (see [i3| for details). We fix eo, io 
and the freeze-out temperature to the same values as for 
a central Au-Au collision at RHIC with Glauber initial 
conditions [10], and p = 3 fm, corresponding roughly to 
the rms values of x and y. Unless otherwise stated, re- 
sults are shown for pions at freeze-out. Corrections due 
to resonance decays [Sy] are not included in this section. 
They are included only in our final predictions in Sec. lIVI 
The second model is a relativistic Boltzmann equation for 
massless particles in 2+1 dimensions (see [ll| for details). 
The only parameter in this calculation is the Knudsen 
number K = X/ R, where the mean free path A and the 
transverse size R are defined as in fll] . p in Eq. ^ is the 
rms width of the energy distribution, while R is defined 
from the rms widths ax and (jy of the particle distribution 
by R~^ = (j~^ + a~^ . For a two dimensional ideal gas of 
massless particles, the particle density n is related to the 
energy density through n ex e^'^ , which gives R = ^p- 
Boltzmann transport theory is less realistic than hydro- 
dynamics for several reasons: 

• the equation of state is that of an ideal gas, while 
the equation of state used in hydrodynamics is 
taken from lattice QCD: it is much softer around 
the transition to the quark-gluon plasma. Al- 
though transport is equivalent to ideal hydrody- 
nanfics when the mean free path goes to zero, our 



results from transport and ideal hydrodynamics dif- 
fer in this limit, because of the different equation 
of state. 

• there is no longitudinal expansion. 

• particles are massless. 

The main advantage of transport theory is that it can 
be used for arbitrary values of the mean free path, while 
hydrodynamics can only be used if the mean free path is 
small. Furthermore, the time evolution of the system can 
be studied and no modeling is required for the freeze-out 
process using transport approach, since one follows all 
elastic collisions until the very last one. 
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calculations [39] . Note that in our hydrodynamic calcula- 
tion, chemical equilibrium is maintained until freeze-out. 
When chemical freeze-out is implemented earlier than ki- 
netic freeze-out, V2/S2 is slightly larger [1^. Fig.[2]shows 
that v^/et, has a magnitude comparable to vil£2i while 
f4/£4 is significantly smaller. Our results for us/es (not 
shown) are even smaller. 
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B. Time dependence 
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FIG. 3: (Color online) Vn/en versus time in transport theory. 
Each curve is the result of a single Monte-Carlo simulation 
with K — 0.025 and £„ — 0.1. The number of particles in the 
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FIG. 2: (Color online) Vn versus En in transport theory and 
ideal hydrodynamics. The Knudsen number in the transport 
calculation is K = 0.025, close to the ideal hydrodynamics 
limit K ^0. 

Fig. [5] displays w„ versus e„ for n = 2, 3, 4 in transport 
theory and ideal hydrodynamics (zero viscosity). The 
values of Vn are smaller in hydrodynamics, which is due 
to the softer equation of state [33] • 

As expected from previous studies of V2 [33| and V3 Q , 
we observe that w„ is linear for small values of £„■ Non 
linearities are stronger for larger values of n, both in 
transport theory and hydrodynamics. A possible inter- 
pretation of these strong nonlinearities is that the the 
contour plot of the initial density is no longer convex if 
e„ > 2/(n^ — 2). The threshold values for n = 3,4 are 
£3 = 1 and 64 — J. If the contour plot is not convex, the 
streamlines (which are orthogonal to equal density con- 
tours) are no longer divergent: shock waves may appear, 
which hinder the development of anisotropies. 

The results presented in the remainder of this section 
are obtained in the linear regime where v„ oc e„. In this 
regime, we find V2IS2 — 0.21, in agreement with other 



In the transport approach, one follows all the trajec- 
tories of the particles, so that Vn is well defined at all 
times, which is not the case in hydrodynamics before 
freeze-out. Fig. [3] displays the results for «„ versus t/p, 
where p is the width of the initial distribution, Eq. (j4]). 
As expected for dimensional reasons [33] , anisotropic flow 
appears for i ~ p. However, w„ appears slightly later for 
larger n. This can be traced to the behavior of Vn at 
early times. The transport results presented in Fig. [3| are 
obtained with a very small value of the Knudsen num- 
ber, K = 0.025, close to the ideal hydrodynamics limit. 
In ideal fluid dynamics, the fluid transverse velocity in- 
creases linearly with t, and w„ involves a n"^ power of 
the fluid velocity, so that u„ scales like t". In transport 
theory, the number of collisions increases like t at early 
times, which gives an extra power of t, and w„ increases 
like i""'""'^ [ll|. In both cases, the behavior of w„ at small 
t is flatter for larger values of n, which is clearly seen in 
Fig. [31 

While elliptic flow keeps increasing with time (it 
slightly decreases at later times, not shown in the fig- 
ure), Vn with n > 3 reaches a maximum and then de- 
creases. The decrease is more pronounced for larger n: 
The mechanism producing Vn is self quenching. 
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FIG. 4: (Color online) Vn/sn versus pt in ideal hydrodynam- 
ics, with En = 0.1. 



Fig. 0] displays the differential anisotropic flow Vn{pt) 
versus the transverse momentum pt for pious in ideal 
hydrodynamics, scaled by the initial eccentricity £„. At 
low Pt, one generally expects w„ to scale like (ft)" for 
massive particles [401 ^. C)ne clearly sees that w„ is much 
flatter at low pt for larger values of n. For larger values 
of Pt , Vn {pt ) is linear in pt . The arguments that explain 
this linear dependence for V2 [l2l can be generalized to 
arbitrary n [41| . The linear behavior at larger pt is also 
clearly seen in Fig.|4l It has already been noted for v^ 9]. 

The value of v^ increases with pt , which explains why 
the ridge and shoulder are more pronounced with a 
high pt trigger ( "hard" ridge) '42| . Though the relative 
strength of vs , is smaller at low pt , it is still comparable 
to V2 , leading to the smaller "soft" ridge [43| . Predictions 
for vslpt) in viscous hydrodynamics for identified parti- 
cles are presented in Sec. IIVI 



D. Viscous damping of w„ 

We study the effect of viscosity first in the transport 
approach, then in viscous hydrodynamics. In trans- 
port, the degree of thermalization is characterized by 
the Knudsen number K . Experimentally, 1/K scales like 
{l/S){dN/dy), where dN/dy is the multiplicity per unit 
rapidit y, a nd S is the overlap area between the colliding 
nuclei f44!|. The dependence of «„ on K can be studied 
by varying the collision system and the centrality of the 
collision 14511 



^ There is no such constraint for massless particles where the 
pt — > limit is singular. Our transport calculations for mass- 
less partons give v„{pt) oc pt at low pt for all n. 
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FIG. 5: (Color online) i;„/e„ versus l/(n^_ft') in transport 
theory. Values of £„ are £2 = ^, £3 = £4 = 0.3 and £5 — 0.1. 
Arrows indicate our extrapolation to K = Q (ideal hydrody- 
namics limit) using Eq. (|6]). 



Transport is equivalent to ideal hydrodynamics in the 
limit X ^^ 0. For small K, observables (such as w„, or 
particle spectra) deviate from the K = ^ limit by correc- 
tions which are linear in K . These are the viscous correc- 
tions: both K and the shear viscosity r\ are proportional 
to the particle mean free path A. Viscous damping is ex- 
pected to scale with the wave number k like k'^ . Here, the 
wavelength of the deformation is 2'iiR/n, hence k ~ n/R. 
Therefore viscous corrections should scale with K and n 
approximately like n^K [4g|. The limit K ^- oo (free 
streaming) is also interesting, since u„ vanishes in this 
limit. For large K, one therefore expects «„ to scale like 
1/K, which is essentially the number of collisions per 
particle [11|. For intermediate values of K {K ^ 1), no 
universal behavior is expected, and observables depend 
on the scattering cross section used in the transport cal- 
culation (dependence on energy and scattering angle). 

Fig. [5] displays the variation of ti„/e„ versus the scaling 
variable l/iji^K) in the transport calculation. Our nu- 
merical results can be fitted by smooth rational functions 
(Fade approximants) [43] for all K: 
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1 + (A„ + Bn)K + C„if 2 + EnK^ - 
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where v^^. An, Bn, C„, D„ and En are fit parameters. 
This formula has the expected behavior in both K ^ 
and K ^)- 00 limits. For n — 2, the lowest-order formula, 
with B2 = C2 = D2 = E2 = 0, gives a good fit ^. For 
n = 3, we obtain a good fit with using the next-to- leading 
order approximant, with D^ = E^ — but free B^, C3. 
For n = 4 or 5, we need all 6 parameters to achieve a 
good fit. Fits are represented as solid lines in Fig. [SJ 
and extrapolations to K — are indicated by arrows. 
As already noted above, the hydrodynamics limits v^^/sz 
and v^2 1^2 are comparable, while v^^ /e^ is smaller by 
roughly a factor of 2. v^^ /e^ is found to be further smaller 
by about a factor 5, with a large theoretical uncertainty. 



For small K, Vn{K) 



(l — AnK): the parameter A„ 



measures the magnitude of the viscous correction. Our fit 
gives A2 = 1.4±0.1 [H, A3 =4.2±0.3, yl4 = 11.0±0.9. 
The error bar on Aq is too large to extract a meaningful 
value. For n — 2,3, 4, we observe An oc n" with a = 2.8± 
0.2, closer to n^ than to the expected n^. The fact that 
viscous corrections are larger for larger n also implies that 
the range of validity of viscous hydrodynamics is smaller 
for Vn with n > 3 than for V2. Even after rescaling K 
by n^, corrections are linear in K only for very small K , 
which is why higher-order Fade approximants are needed. 
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FIG. 6: (Color online) Vn/e„ versus rj/s in hydrodynamics. 
The initial and freeze-out temperature are Ti = 340 MeV and 
Tf = 140 MeV, respectively. 



The magnitude of viscous effects can be seen more di- 
rectly by varying the shear viscosity 77 in viscous hydro- 
dynamics [48|. For each value of n, we have performed 
three calculations with 77 ~ (ideal hydrodynamics), 
rj/s = 0.08 ~ l/47r Q, and t]/s = 0.16, where s is 
the entropy density. The result is presented in Fig. [HI 
The variation of w„ with rj is found to be linear for all n 
for this range of viscosities, which is a hint that viscous 
hydrodynamics (which addresses first-order deviations to 
local equilibrium) is a reasonable description. Interest- 
ingly, the lines are almost parallel, which means that the 
absolute viscous correction to w„/e„ depends little on n. 
However, since Vn/Sn is smaller for larger n, the rela- 
tive viscous correction is larger for larger n. From the 
transport calculation, we expect that the relative viscous 
correction is 3 times larger for W3 than for V2 , and 8 times 
larger for U4 than for V2 ■ The increase in Fig |6] is more 
modest. Note that we keep the freeze-out temperature 
constant for all values of 77/s. Strictly speaking, this is 
inconsistent. Freeze-out is defined as the point where 
viscous corrections become so large that hydrodynam- 
ics breaks down: when the viscosity goes to zero, so does 
the freeze-out temperature [12l | . By varying only rj/s and 
keeping Tf constant, we only capture part of the viscous 



correction ^. Since triangular flow, like elliptic flow, de- 
velops at early times, W3 is sensitive to the value of 77/s 
at the high-density phase of the collision. 



IV. PREDICTIONS FOR V3 AT RHIC AND LHC 

A. Triangularity fluctuations 

We now give realistic predictions for V3 at RHIC and 
LHC. The transport calculations in Ref. ;9] show that 
even with lumpy initial conditions, V3 in a given event 
scales like the triangularity £3. We define £„ as in [9|: 



Cr},e 



in^j 



_ J e{x,y)r'^e^"'^dxdy 
J e{x, y)r'^dxdy 



(7) 



are 



where e{x, y) is the initial energy density and (r, 
the usual polar coordinates, x = r cos (j), y = r sin </>. 

Following the discussion in Sec. [Hi experiments mea- 
sure the average value of (wn)^, so that 



V7'' = V(KP)- 
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Assuming u„ — Ke„ in each event, the measured w„ scales 
like the root mean square e„ defined by 



Vii^nF) 



(9) 



We compute e™*' using two different models. The first 
model is the PHOBOS Monte-Carlo Glauber model ;5C|, 
where it is assumed that the initial energy is distributed 
in the transverse plane in the same way as nucleons 
within colliding nuclei. We modify the initial model 
slightly [33| by giving each nucleon a weight w = 1 — 
X + a;iVcou, where iVcou is the number of binary collisions 
of the nucleon. We take x = 0.145 at RHIC and x = 0.18 
at LHC [m. The second model is the Monte-Carlo KLN 
model of Drescher and Nara [52| , which is the only model 
incorporating both saturation physics and fluctuations. 
Both of these models yield event-by-event eccentricity 
fluctuations, which are consistent with measured elliptic 
flow fluctuations [3J| . We loosely refer to the two models 
as Glauber and Color Glass Condensate (CGC). 

Fig. [7] displays er™" as a function of the number of 
participants, e™** is larger than e^s for non-central 
collisions, which is due to the almond shape of the over- 
lap area. The eccentricity is somewhat larger with CGC 
than Glauber [5J]. e™ is very close to e™. Both vary 
with A'^Part essentially like (A^Part)"^^^, as generally ex- 
pected for statistical fluctuations [55|. Unlike e™, they 
are slightly smaller for CGC than for Glauber. Since 
the only source of fluctuations that is considered in both 



^ We have checked that v-^/ez is larger with a lower freeze-out tem- 
perature Tf = 100 MeV. In particular, we find v^/ez > i'2/e2, 
in agreement with the transport calculation. 
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FIG. 7: (Color online) Root mean square eccentricities e™^ 
for n = 2, 3, 4, 5 for Au-Au collisions at 200 GeV per nucleon, 
versus the number of participant nucleons A^part- Apart is 
used as a measure of the centrality in nucleus-nucleus colli- 
sions: it is largest for central collisions, with zero impact pa- 
rameter [53]. Thick lines: Monte-Carlo Glauber model [5CI |: 
Thin lines: Monte-Carlo KLN model fsl. 



models is the position of the nucleons in the colliding 
nuclei, this difference may be due to the technical im- 
plementation of the Monte-Carlo KLN model. Finally, 
e™** is slightly larger than odd harmonics for peripheral 
collisions because the almond shape induces a nonzero £4 
as a second order effect. Fig. [7] only displays results for 
Au-Au collisions at RHIC. Results for Pb-Pb collisions at 
the LHC are similar, except for a different range in A'part, 
and a somewhat larger difference between Glauber and 
CGC for £3. 



B. Method for obtaining va in hydrodynamics 

In order to make predictions for W3, we start from a 
smooth initial energy profile e(r, 0), possessing the usual 
symmetries </> — > — and </> — S- + tt. We then put by 
hand a cos(3(/)) deformation through the transformation, 
inspired by Eq. Q, 



eir 



-> e 



1 



;cos(3(0- Va)),^ 



(10) 



where £3 is the magnitude of the deformation, and f/'a 
the flat axis of the triangle. We choose £3 = e™^. The 
choice of "03 is arbitrary. The initial profile has a nonzero 
eccentricity for noncentral collisions, due to the almond 
shape of the overlap area. Through Eq. (TTO]) . we add a 
triangular deformation to an ellipse. Since the original 
profile has (j) ^ (j) + tt symmetry, -^g is equivalent to 
-03 + f. Furthermore, -03 is equivalent to —-03 due to 
— >■ — symmetry. Therefore, one need only vary ■03 
between and |. We choose the values 0, f^ and ^. 

We then compute £3 and 03 defined by Eq. ([T]). With 
the gaussian profile (JH), the input and output values are 



identical: £3 = £3, 03 = -03. Our predictions use two 
sets of profiles which both describe RHIC data well |l3 : 
optical Glauber and (fKLN) CGC. With both profiles, e'^ 
differs from £3 by a few percent. 03 is essentially identi- 
cal to 03 , which means that the elliptic deformation does 
not interfere with the triangular deformation. According 
to the previous discussion, we should tune £3 in such a 
way that £3 = £™^ in order to make predictions for D3. 
It is however easier to use the proportionality between 
V3 and £3: one can then do the calculation for an arbi- 
trary £3, and rescale the final results by £™/£3. We use 
£™'^ from the Monte-Carlo Glauber model with Glauber 
initial conditions, and from the Monte-Carlo KLN model 
with CGC initial conditions. 

There is some arbitrariness in the definition of the tri- 
angularity £3: one could for instance replace r^ by r'^ in 
Eq. ([7]) [5g|. With this replacement, both £3 and £™ 
(from the Monte-Carlo calculations) increase, but the 
ratio £3™''/£3 — and therefore also our predicted V3 — 
changes little (less than 7% for all centralities and both 
sets of initial conditions)"^. 

Finally, we compute ^3 in viscous hydrodynamics. It 
has been shown that RHIC data are fit equally well with 
Glauber initial conditions and rj/s = 0.08 or with CGC 
initial conditions and rj/s — 0.16 [l3|. The larger eccen- 
tricity of CGC (which should produce more elliptic fiow) 
is compensated by the larger viscosity (larger damping 
and less fiow) , so that the final values of V2 are very simi- 
lar. For LHC energies, details are as in Ref. [53| (with V3 
calculated from a Cooper- Frye freeze-out prescription). 
In all cases, v^, is found to be independent of the orien- 
tation of the triangle tp'^. In the case of Glauber initial 
conditions, we perform calculations of V3 with and with- 
out resonance decays at freeze out [3y| . Resonance decays 
roughly amount to multiplying V3 by 0.75 at RHIC, and 
by 0.83 at LHC. Our CGC results are computed without 
resonance decays, and multiplied by the same factor at 
the end of the calculation. 



C. Results and comparison with data 

Results are displayed in Fig. [5] for both sets of initial 
conditions. CGC initial conditions have both a smaller 
triangularity, and a larger viscosity, so that they pre- 
dict a much smaller V3. The change in viscosity explains 
roughly 70% of the difference between CGC and Glauber 
at RHIC, and about half at LHC. The centrality depen- 
dence is much flatter in Fig. [8] than in Fig. [Tj The de- 
crease of £3™'' with increasing A^part is compensated by 
the increase of the system size and lifetime, which leads 



^ If one replaces r^ by r*^ in Eq. l[7j, £„ scales with k like k + 2 for 
a smooth, symmetric density profile e(r) deformed according to 
Eq. moil . Therefore, £3 is larger by ^ if defined with a factor r^ 
instead of r^ . 
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FIG. 8: (Color online) Average W3 of pions as a function of the 
number of participants for Au-Au collisions at 200 GeV per 
nucleon (RHIC) and Pb-Pb collisions at 5.5 TeV per nucleon 
(LHC). Hydrodynamic predictions are for Glauber initial con- 
ditions with Tj/s = 0.08, and CGC initial conditions with 77/s 
= 0.16, which best fit 112 data at RHIC fTot]. 



to a smaller effective Knudsen number K or, equivalently, 
a smaller viscous correction. We predict values of v^ sig- 
nificantly larger at LHC than at RHIC. This is because 
viscous damping is less iniportant due to the larger life- 
time of the fluid at LHC [57| . 

Although experimental data for triangular flow are not 
yet available, both V2 and V3 can be extracted from 
the measured two-particle azimuthal correlation using 
Eq. ^ [9]. Figs. [^ and [TUl display a comparison between 
experimental data from STAR ^] and our hydrodynamic 
calculations. The STAR data is obtained from correla- 
tions between particles at midrapidity (|?7| < 1) and in- 
termediate transverse momentum (0.8 < pt < 4.0 GeV). 
The correlation results have been projected at 1.2 < 
A77 < 1.9 to reduce sensitivity to nonflow effects. 

We first discuss our results for V2- As explained above, 
our hydrodynamic model has smooth initial conditions, 
and does not include the effect of eccentricity fluctuations 
for V2- Since V2 oc £2 to a good approximation, we have 
rescaled our result for V2 by the rms £2 from Fig. [7] (again 
using the Monte-Carlo Glauber for the Glauber initial 
conditions and the Monte-Carlo KLN for CGC). This 
rescaling significantly improves the agreement with data, 
compared to ,10], for the most central bin. As shown 
in Fig. |9l the agreement between theory and data is ex- 
cellent with both sets of initial conditions. The smaller 
viscosity associated with Glauber initial conditions re- 
sults in a somewhat steeper centrality dependence than 
for CGC initial conditions. 

Results for V3 are shown in Fig. [TUl The larger magni- 
tude, compared to Fig. [51 is due to the low pt cutoff. The 
cutoff also enhances the effect of viscosity, resulting in a 
larger difference between Glauber and CGC. With a low 
Pt cutoff, the viscous correction is mostly due to the dis- 
tortion of the momentum distribution at freeze-out ISSlI . 



FIG. 9: (Color online) 112 for charged particles with 0.8 < 
Pt < 4 GeV/c extracted from STAR charge-independent cor- 
relation data |4|, and predictions from viscous hydrodynam- 
ics \1^ with Glauber initial conditions and 77/s = 0.08, or 
CGC initial conditions with rj/s — 0.16. Theoretical calcula- 
tions are for pions with the same pt cut as data, and scaled 
by the rms eccentricity from the corresponding Monte-Carlo 
model. See text for details. 
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FIG. 10: (Color online) Same as Fig.[9]for V3. 



The momentum dependence of this distortion is strongly 
model-dependent [59|. The present calculation uses the 
standard quadratic ansatz, which may overestimate the 
viscous correction at large pt [60|. The magnitude and 
the centrality dependence of V3 observed by STAR are 
rather well reproduced by our calculation with Glauber 
initial conditions, except for peripheral collisions where 
hydrodynamics is not expected to be valid. 

Fig. [TT] displays our predictions for ^3(^4) of identified 
particles at RHIC. As anticipated in Ref. J41|], the well- 
known mass ordering of elliptic flow |i6l| is also expected 
for V3. At high pt, a strong viscous suppression is ob- 
served. As explained above, the pt dependence of the 
viscous correction is model dependent, and it is likely 
that the quadratic ansatz used here overestimates the 
viscous corrections at large pt [60|. Note that effects of 
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FIG. 11: (Color online) Differential triangular flow for iden- 
tified particles in central (0 — 5%) Au-Au collisions at RHIC. 



resonance decays are not included in Fig. [TT] Resonance 
decays only change the results slightly in the low-pt re- 
gion. 

Finally, we have also computed v^ along the same lines 
as U3. The driving force for v^ is the rms £5, which is 
very close to £3 (see Fig. [7]). However, the hydrodynamic 
response is much smaller, and viscous damping is also 
much larger as discussed in Sec. IIIII We find that the 
average integrated vq is smaller than W3 by at least a 
factor of 10. Results for differential v^ are presented in 
Fig. 1121 vq varies more strongly with pt than V2 and V3, 
and becomes as large as 0.02 at pt = 1.5 GeV/c if the 
i]/s is as small as 0.08. For larger viscosity, however, v^ 
may be too small to measure even with a high pt trigger. 




FIG. 12: (Color online) Differential pentagonal flow for iden- 
tified particles in central (0 — 5%) Au-Au collisions at RHIC. 



V. CONCLUSIONS 

We have presented a systematic study of triangular 
flow in ideal and viscous hydrodynamics, and transport 
theory. Triangular flow is driven by the average event- 
by-event triangularity in the transverse distribution of 
nucleons, in the same way as elliptic flow is driven by the 
initial eccentricity of this distribution. The physics of W3 
is in many respects similar to the physics of V2 ■ In ideal 
hydrodynamics, the response to the initial deformation is 
almost identical in both harmonics: v^/e^ ~ W2/£2 — 0.2. 
For quadrangular flow, Vi/en is smaller, typically by a 
factor 2. For pentagonal flow, us/es is so small that U5 
is unlikely to be measurable, even though £3 and £5 are 
almost equal. W3 develops slightly more slowly than W2, 
though over comparable time scales. The dependence on 
transverse momentum pt is similar for v^ and V2, but 
vz/v2 increases with pt- Hydrodynamics predicts a simi- 
lar mass ordering for v^(pt) and V2(pt)'- i'3 at fixed pt is 
smaller for more massive particles. These results can be 
checked experimentally by a differential measurement of 
triangular flow. 

We have also made predictions for triangular flow, 1)3, 
at RHIC and LHC, using viscous hydrodynamics. Using 
as input the triangularity from a standard Monte-Carlo 
Glauber model, and a viscosity 77/s = 0.08, we repro- 
duce both the magnitude (within 20%) and the central- 
ity dependence of W3 extracted from STAR correlation 
measurements, without any adjustable parameter. Our 
results support the hypothesis made in Ref. @ that trian- 
gular flow explains most of the ridge and shoulder struc- 
tures observed in the two-particle azimuthal correlation. 

Triangular flow is a sensitive probe of viscosity. Vis- 
cous effects drive the energy and centrality dependence of 
f3. More central collisions have less fluctuations, hence 
smaller triangularity. This decrease is to a large extent 
compensated by the increase in the system size and life- 
time, resulting in a very slow decrease of W3 with cen- 
trality (except for peripheral collisions where viscous hy- 
drodynamics is unlikely to be valid). Comparison with 
existing data favors a low value of 77/s. At LHC, smaller 
viscous corrections are expected due to the increased life- 
time of the fluid: we predict that v^ should be larger than 
at RHIC, typically by a factor |. 

The absolute value of W3 scales linearly with the aver- 
age initial initial triangularity. We have used two models 
of initial geometry which incorporate fluctuations, the 
Monte-Carlo Glauber model and the Monte-Carlo KLN 
model. The underlying source of fluctuations is the same 
in both of these models. More work is needed to constrain 
initial fluctuations on the theoretical side. More work is 
also needed to incorporate these fluctuations more read- 
ily into hydrodynamic calculations. Although triangular 
flow is expected to be created by lumpy initial conditions, 
our predictions are based on smooth initial conditions, 
in the same spirit as the study of transverse momentum 
fluctuations of Ref. [22|. The underlying assumption is 
that vs/es is the same for lumpy initial conditions and for 



smooth initial conditions. The vahdity of this assump- 
tion should eventually be checked. 

Triangular flow is a new observable which should be 
used to constrain models of heavy-ion collisions, along 
with elliptic flow. Elliptic flow depends on initial ec- 
centricity, fluctuations, and viscosity, which are poorly 
constrained theoretically. Triangular flow solely depends 
on fluctuations and viscosity, with a stronger sensitivity 
to viscosity than V2 ■ Two different sets of initial condi- 
tions, which fit V2 data equally well, give very different 
results for v^. Experiments could measure W3 as a func- 
tion of transverse momentum, system size and centrality. 
As shown in this paper, theoretical predictions for the de- 
pendence of U3 on these parameters are very specific. If 
experiments confirm our predictions, simultaneous anal- 
yses of V2 and W3 can be used to improve our understand- 
ing of the initial geometry of heavy-ion collisions, and pin 



down the viscosity of hot QCD. 
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